ON SCHRODINGER EQUATION WITH TIME-DEPENDENT QUADRATIC 

HAMILTONIAN IN R d 



ERWIN SUAZO 

Abstract. We study solutions to the Cauchy problem for the linear and nonlinear Schrodinger 
equation with a quadratic Hamiltonian depending on time. For the linear case the evolution operator 
can be expressed as an integral operator with the explicit formula for the kernel. As a consequence, 
conditions for local and global in time Strichartz estimates can be established. For the nonlinear 
. case we show local well-posedness. As a particular case we obtain well-posedness for the damped 

harmonic nonlinear Schrodinger equation. 



1. INTRODUCTION 

In this paper we first study the time-dependent linear Schrodinger equation 



i^ = H(t)^ 1>(x,0) = <p (1.1) 
J> ' 

^ ■ with the quadratic Hamiltonian 

<N ; H(tW = + £ - + %(*)|^ - (tej^ +*)), (1-2) 

psj . i=l ^ 3 \ 3 / / 

where bj, fj,gj, Cj G C 1 (bj, fj, Qj could be piecewise continuous functions) and <p G S(M. d ) (S{R d ) is 
the Schwartz space) to simplify the discussion. We derive its evolution operator given by an explicit 
J> ! formula in the form 

X'- iP(x,t) = U H {t,0)cp{x) = / G H (x,y,t) <p(y) dy (1.3) 

by constructing the fundamental solution (FS) associated to (jl.ip - (jl.2p . see Lemma 1; the FS is a 
solution of (11. ip with the initial data Gn(x,y,0) = 8(x — y). Next we study properties of Uh and 
give conditions to obtain local and global in time Strichartz estimates for ip G L 2 (M. d ). We prove that 
the nonlinear version of f ll.ip -f fl~2]) with algebraic nonlinearity is locally wellposed in L 2 x (R d ) in the 
subcritical sense, see Section 4. Finally, we introduce the damped harmonic nonlinear Schrodinger 
equation: 

.du un ( d 2 u , \ .A / du \ , , .„ , , s 
i— = — -— + x 2 u + i- 2x— + u + h\u\*- x u (1.4) 



ft 2 \ & 2 J 2 V dx 

and also prove well-posedness. A systematic study of the blow up and scattering results of fll.ip - (ll.2p 
are presented in [52"] . 



1991 Mathematics Subject Classification. Primary 35A08, 35Q40, 35Q55. Secondary 81Q05. 

Key words and phrases. The Cauchy initial value problem, Schrodinger equation (mathematical properties) , Qua- 
dratic Hamiltonian, Fundamental solution, Propagator, Green function, Mehler's formula, Strichartz estimates, Non- 
linear Schrodinger equation. 

1 



2 



ERWIN SUAZO 



The study of time dependent quadratic Hamiltonians is quite complex and has brought a great 
deal of attention see [21] , [25] and references therein, see also [22] for a nice discussion on this type 
of Hamiltonians consider "folk wisdow" . The study of methods (i.e. propagator method, dynamical 
invariant and second quatization methods) used to find the exact propagators for Schrodinger 
equations has also been studied by several authors [55]. The equation fll.ip can be solved, at least 
formally, using a time evolution operator, U (t, t ) given by 

U (t, t ) = jy~ l ) k [ ^i- I ^ dt k H{t x )...H{t k ) = T (exp (-L H (f) dt'^j ) , (1.5) 

where T is the time ordering operator which orders operators with larger times to the left, and of 
course this expression might diverge. 

The fundamental solution Gh for the equation fll.ip - fll.2l) includes the following examples with 
the explicit expressions: 



Table I. Some exactly solvable quadratic Hamiltonians (We assume E, k constants). 



Hamiltonian H(t) 


Fundamental Solution (Propagator) 


Free Particle 


G (x,y,t) - f -l—S x -y\ 2 l 2t 

uv ' ' v 2~i smt 


Constant Electric Field 


xexp^f^t f 4 V) 


Isotropic Oscillator 
H 2 {t)^ = -\ d ^ + \x 2 ^ 


G 2(x,y,t) 

x ex P (^sit) (( x2 + y 2 ) cost 2x y)) 


Repulsive harmonic potential 

W)l> = -h& - l*V 


Gs(x,y,t) — . 1 . ; 

xex P ( z 4 Si nh j(a; 2 + l/ 2 )cosht 2xy)) 


Anisotropic Oscillator 


GJx,y,t) - - w . - 

X eX P ( 2 4sinM) ^ + C ° S Ut 2Xy }) 




Modified MC-SS Oscillator 
H e (t)il) = - cos 2 tft + sin 2 tx 2 ip 
-i^ (2s£ - 1) * 


Gq (x, y, t) — , 1 

\/ 27ri(cos t sinh t+sin t cosh t) 
I [x 2 — j/ 2 ) sint s\nht+2xy— (ie 2 +j/ 2 ) costcosht \ 


exp 1 2i (cos t sinh t+sint cosh i) ) 


Damped Harmonic Oscillator 

H 7 (t)i/> = -<f%£ + fx 2 ^ 

+il (2x1 + 1) i> 


G 7 (x,y,t) = J . -. t 

X eX P f 2^no,t CO** + ^) COSWt 2 ^) y 

x exp (x 2 - y 2 )) , c = v/c 2 - A 2 : 


I 

> 


Analog of Heat Equation 
with Linear Drift 

H 8 m = -^-ikx§, k>o 


^8{x, y ,t) ^ iainh(kt) c ^y 4smh( fci) y 



The generality of equation (ll.l[) -( TTT2|) includes several examples of interest not only for linear [18J, 
[55] , but for nonlinear Schrodinger equations (NLS) [5], [6]; for a general review in NLS see [Hj, [56J. 
The explicit solution formula for fll.ip ~ fll.2p . see (I2.10p . allows us to treat nonlinear versions of these 
cases with the advantage of allowing time dependence factors so that questions such as local and 
global existence, finite time blow up and scattering can be solved using similar methods to those 
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used in the case of the nonlinear Schrodinger equation without potential. The results presented here 
are of interest to study the behavior in time of singularities of solutions for Schrodinger equations 
in the same direction as in |61j. Also the fact that in quantum mechanics it is rare to find an exact 
solution to nonstationary problems, see [43], [12], makes our explicit solution useful for testing 
numerical methods to solve time-dependent Schrodinger equation. One of the original motivations 
of this paper is to introduce a generalization of the practical formula fjl.19p - fjl.20p . for applications 
see [5J, [61] and [33] . 

1.1. Quadratic Quantum Hamiltonians. The expert would recognize fjl.2j) as a quantum me- 
chanical self-adjoint Hamiltonian, which is a quadratic polynomial in x and p = —id/dx with time 
dependent coefficients 



H(t) = a(t)p 2 + C M(p. x + x .p) + h M x * - g (t)p - f(t)x + C(t). (1.6) 

As pointed out in [23] one can assume ((t) = since it causes a trivial phase factor in the propagator. 
We also have assumed in (II .ip a(t) = 1/2, thinking we can disregard it after a substitution, however 
there are important cases that we will be missing, for example the Caldirola-Kanai Hamiltonian [I] , 
|29j . In [12] the authors did not know about this case and called it "third model," it was one the 
models of the damped harmonic oscillator with explicit propagator considered in that publication. 
The Caldirola-Kanai Hamiltonian was introduced more than 60 years ago [I], [29] : 

HcK(t)ip — H x 2 ib. (Caldirola-Kanai Hamiltonian) 

2 ox 2 2 

its fundamental solution is given by 



where 



G C K(x,y,t) = J UeXt . e^+mzv^)^ (L7) 
Z7TtUo sin out 



uj cos out — A sin out 9i , 

a(t) = e 2Xt , (1.8) 

2cuq sin out 

0{t) = ^-e Ai , (1.9) 

ouq sin out 

ou cos out + A sin cut 

7(t) = . 1.10 

,w 2ou smout v ; 

This model has been studied for several authors, and a detailed bibliography can be found in [55] . 

Another example with a(t) not constant that can be solved using the same ideas presented in 
this work is 

H 10 (t)^ = -a(t)^ + b -^-x 2 ^ (1.11) 

= (^ 2 cos(ta)- 7 sin(ta) tanh( 7 t)) 
cosh(7t)(cos(7i) cosh(7t) — 27) 

ou 2 - 



m - - Wy n = V^-^, (1.13) 
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and its fundamental solution is given by 



Gio(x,y,t) = Ho^cosh7t t(o( ^ 2+/3ft)w(%2); 



with 



cosh(7t) (m fi cosh(7t) cos(fii) - 7) 
= ' (L15) 

= _ m ncosM (1 16) 

2 sm(JZr) 

Suslov et. al in [13J study the quantum integrals of motion between others of these last two models. 
1.2. A General Formula. If we consider the linear Schrodinger equation on R n (ll.ip - (ll.2l) with 

H v (t)1> = ~A1> + V(x)i/>, (1.18) 

where V(x) = Y^j=\ 0>j~2~ x j + n — 1> w i > 0' e { — 1; 0' 1}> then the solution is given by 

the following formula [5] 

iK*,t) = Uv(t)f:= /— - — — / ^y^f{y)dy (1.19) 
y/2tngj{tj Jr™ 

where 

'-/'.;(') •••'•.///./ 




£15^2^1 ; cosh(a;jt)|, if <5j = — 1, 

{t, 1}, if = 0, (1.20) 
^,cos(c^)), if ^ = +1. 



This formula includes the free particle propagator introduced by Ehrenfest [TF] that corresponds 
to the case 8j = for all j in (I1.18p . The propagator for the simple harmonic oscillator is obtained 
from (11.191) . by choosing 8j = 1 for all j in (I1.18P (a consequence of Mehler's formula for Hermite 
polynomials [18]). Finally, the propagator for the isotropic harmonic potential is obtained by 
choosing 8j — — 1 for all j in (11.181) . 

The fact that we want to emphasize in formulas (I1.19l) - (ll.20l) is that it combines the fundamental 
solutions of the free particle and simple and isotropic harmonic oscillators in the one dimensional 
case (using tensor product) to construct the explicit solutions from a variety of Schrodinger equations 
in several dimensions. One motivation comes from the Schrodinger equation (ll.ll) - (ll.2p with the 
Hamiltonian 

Hh(t)i> = -g AV> + 1/2 {-ujix\ + u 2 x 2 2 ) </>, (1.21) 

(compare with H 3 and H4 from Table I) whose evolution operator satisfies global in time Strichartz 
estimates, see [5]. 
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However, the explicit formula for the evolution operator corresponding to fll.ip - fll.2l) allows us to 
consider a wider range of operators to solve explicitly a variety of Schrodinger equations in several 
dimensions. 

1.3. Cases from Table I Not Included in Formula ( ll.l9j) -( Tl.20j) . The case of the particle in 
a constant external field having the Hamiltonian H\(t) case was studied in detail by [15], [JJ, [3], 
[18] . [25] . |4T] and [45] . The forced harmonic oscillator is obtained in f ll.ip - fll.2l) by choosing 

where S(t) is a complex valued function of time t and the symbol * denotes complex conjugation. 
It corresponds to the case which is of interest in many advanced problems, examples including 
polyatomic molecules in varying external fields, crystals through which an electron is passing and 
exciting the oscillator modes, and other interactions of the modes with external fields. An example 
of this type of Hamiltonian is H 6 (t), introduced in [TO] . 

The damped oscillations have been analyzed to a great extent in classical mechanics; see, for 
example, [2] and [31]. In [T2] the damped harmonic Schrodinger equation with self-adjoint Hamil- 
tonian Hr(t) was considered. The case H$(t) is taken from the analogous case of the heat equation 
with linear drift, see [38] for a proof of this fundamental solution using Lie theory approach. 

All examples in Table I are explicit solutions because one can solve the characteristic equation 
(12. 2p and evaluate the integrals (I2.5p - fl2.8l) explicitly. Solving (12. 2p could require elaborated tech- 
niques, see for example [35]. However, even without solving (I2.2p explicitly important conclusions 
can be drawn. For example, see [7J, where 

Hy(t)ij = -^AtP + V(x,t)^ (1.22) 

with V{x, t)il) = J2 d j=1 hj ^x 2 ^ was considered and finite time blow up results are obtained. Another 
interesting example is the parametric forced harmonic oscillator with a time- dependent frequency 
u(t), see 02] 

! 777. 

H (t)iP = + j0j(t) 2 x 2 iP - f(t)xip. (1.23) 

This case is relevant in the context of charged particle traps [21] . 

1.4. Organization of the paper. In Section 2 we derive the explicit formula for the fundamental 
solution of (11. ip with Hamiltonian (11. 2p (Lemma 1). We also present explicitly the time evolution 
operator and relevant properties (Corollary 1), conditions on the uniqueness of the solution for (11.11) - 
(11. 2p . and continuous dependence on the initial data and smoothness of the solution (Theorem 1). 
In Section 3 we discuss Strichartz type estimates for Uh (Lemma 5 and Theorem 2). Finally, in 
Section 4 we discuss local well-posedness in L 2 .(M. d ) for the nonlinear case in the subcritical sense 
(Proposition 1). 

2. LINEAR CASE 



In this section we wish to prove a generalization of the formula fll.19p - fll.20l) in M d : 
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Lemma 1. (Fundamental Solution) 1. Let tp e S(M. d ). The Cauchu initial value problem I[1J\) - U.2\) 
has the following fundamental solution 

/ d 1 \ d / 2 

G H (x,y,t)= I Y[ J e *(£M*)z 2 +/M i )^'W+7 J (%^ (2.1) 

)Ltj satisfies 

fjfj + Acjj (t) n s = 0, (2.2) 

witt Oj (t) = bj(t)/2 — Cj(t)/4 — c'j(t)/4:, which must be solved subject to /^(O) = 0, /Xj-(O) = 1. 
Furthermore, ctj (t) , /3j (t) , jj (t) , 8j (t) , Sj (t) , Kj (t) are differentiable in time t only and are given 
explicitly by 

ai (i) = if£M _ sW (2 . 3) 

1 w 2 ^ (t) 2 ' V ; 



1 



ft ( £ ) = " — . ( 2 - 4 ) 



2^(0 ^(t) Jo (^.(r)) 2 2 1 ; 



3 

t 



6j (t) = -L, / (/, (r) - c, (r) # (r)) (r) + % (r) /jj (r) dr, (2.6) 
A*j l r J Jo 

^ (^.(r)) 2 

i 

Jo V ( r ) 

*i (*) = MU? (*) " 2 T fa (r) 5 S (r)) 2 rfr (2.8) 

Vi( ^ (T) (/i(r)-c,(r)^(r)) rfr 
o ( r ) 



with 



5 3 (0) = # (0) , £j (0) = (0) , Kj (0) = 0. (2.9) 
2. Convergence to the initial data: 



lim / G H (x,y,t)ip(y,t)dy = ip(x,0). 

Thus, the fundamental solution (propagator) is explicitly given by (I2.14p in terms of the charac- 
teristic function (Q with (TO) - fl2^B|) . 

Remark 1. The conditions A2.9\) , which are justified since we are looking for the following asymp- 
totic formula, see [52], holds: 



->■ -== exp % 3 3 I exp (i gj (0) (a^ - Vj )) 



A /27ri/7~(i) V 2t 
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Theorem 1. 1. Let <p G S(R d ), then the Cauchy initial value problem U.1\) - /[T72)) has the following 
unitary evolution operator: 

( d i \ d ' 2 r 

U H {t)ip=\Y\ / e^ ^?^'^*^ (2 10) 

2. Iffpe S(R d ), then U H (t)<p G S(R d ). 

3. Un(t,s) = UH(t)U^ 1 (s) and by duality Uh (t, s) can be extended to S'(R d ). Furthermore, 
U H (-)<peC(R, S'(R d )) for every <p G S'(R d ). 

If if) satisfies $l.l\) - (TJfy and it is smooth, then: 

4- The following estimates hold: 

||^ff(XMU 2 (R d ) = IMU 2 (R d )> (2-H) 

( d i y /2 

\\U H (t, s) ^ < [U ^ (t) ^ {s) (5) _ Tj . (t)) J IMIxH*,- ^ 12 ) 

5. Uniqueness and continuous dependence on the initial data in L 2 x (R d ) holds. 

Remark 2. The function fij will characterize the singularities [H], [SU], [EE] &n<i [58] . More 
complicated cases may include special functions, like Bessel, hypergeometric or elliptic functions. 
For a nice conection between the characteristic equation \2. 2\) and Ehrenfest Theorems see 



Corollary 1. The evolution operator associated to ( fi. ij) -( TOj) satisfies the following properties: 

1. U H (t,s) = U H (t)U H \s). 

2. U H {t,t) = Id. 

3. The map (t,s) — >■ Unites) is strongly continuous. 
I U H (t,T)U H (r,s) = U H {t,s). 

2.1. Proof of the Lemma 1. We follow [TT] where the fundamental solution is constructed for 
a more general case in d = 1 dimension. The Lemma follows from the construction in the Id case 
and then using tensor product to construct the <i-dimensional fundamental solution. We recall how 
to construct the fundamental solution in the one dimensional case for 

**' - 1 **' + - fM + <*M^ ^ ■ (2.13) 



9* 2 c^ 2 2 jrj Jn ' jrj ' w 2 V 

The fundamental solution is found by using the ansatz 

V>j = Aj-e*^ = Aj- (t) e^^'^ (2.14) 

with 

A ' =A ' {t)= vmm (2 ' 15) 
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and 

Sj = Sj (x, y, t) = atj (t) x) + (3j (t) xjyj + ^ (t) y) + Sj (t) x 3 - + e 3 - (t) yj + k 3 - (t) , (2.16) 

where <x,- (t) , j5j (t) , jj (t) , 8j (t) , e 3 - (t) , and Kj (t) are differentiable real- valued functions of time t 
only. Indeed, 

— 1 = i^r 1 ) - hx 2 , + fjXj + (o 7 - - CaXa) — - 1 (2.17) 

dt 2 V dx J J j; dxj v 7 

by choosing 

// 19 2 5 7 - c,- , , c,-(t) 

Equating the coefficients of all admissible powers of x™y™ with < m + n < 2 gives the following 
system of ordinary differential equations: 



da ■ 

_1 + 6 . ( t ) + 2cj (t) a 3 + 2a) = 0, (2.19) 

-jZ + (c j (t) + 2a j (t))l3 j = 0, (2.20) 

l + ^r = °- < 2 - 21 > 

+ ( Cj (t) + 2a 3 (t)) (5, = £ (t) + 2a, (t) % (t) , (2.22) 

^ = (9i(t)-8At))ft(t)> (2-23) 

^=*(W*)-^, (2-24) 

where the first equation is the Riccati nonlinear differential equation. Substituting (I2.18P into (I2.19P 
results in the second order linear equation 

fi'j + Aaj (t) fij = (2.25) 

with 



b 3 (t) cf(t) 



= ^-^r-j. ( 2 - 26 ) 

which must be solved subject to the initial data 

H S (0) = 0, /4 (0) = 1. (2.27) 

We shall refer to equation ( 12.251) as the characteristic equation and its solution fij (t) , subject 
to (I2.27p . as the characteristic function. Using integration by parts we can solve (12.1 91) - 02. 241) 
obtaining (F23D-0- 

Part 2 is a consequence of the two following results in the one dimensional case proven in [52J. 
Lemma 2. Let G be defined by Ii2.1\) . There exists a complex-valued function K satisfying 

G(x,z,t)K(z,y,Q) X (y)dzdy= [ [ G{x, z,t)K(z,y,0) X (y)dydz (2.28) 
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and 

[ G{x,y,t)il>(y,t)dy= [ K(x,y,t) X (y)dy. (2.29) 

Lemma 3. If we consider the intial data i/j(x, 0) such that if)(x, 0) = J Rd K(x, y, 0)x{y)dy for some 
X e L 1 ^), then 



lim / G(x,y,t)ip(y,t)dy = ip(x,0). 

Remark 3. For the case of the free particle propagator K(x, y, 0) = e tx ' y and x{v) i> s the Fourier 
transform ofip(x, 0). 

2.2. Proof of Theorem. 1. The propagator Unit) for the equation (I1.10 - (ll.2p can be written as 

U H (t)<p = AtBt8{Ct<p), (2.30) 

where A t (x) = e ^(t)^+s(t>+<t)) ^ = ^{t^+e^x) ^ B t w(x) = (2m fi(t))~ * w {-/3{t)x/ 2ir) and 

$ is the Fourier transform. Since the Fourier transform is an isomorphism on the Schwartz space 
we have the operator Unit) is an isomorphism on the Schwartz space. 

2. We claim that 



^ (x, 0) = U H X (t) if>(x,t)= H [x, y, t) V [y, t) dy, (2.31) 

J — oo 

where 
such that 

U H (t) Uh 1 (t) = Uh 1 (t) U H (t) = I = id. (2.33) 
First we observe that the following orthogonality relations of the kernels hold: 

„ d 

/ G(x } y } t)H(y } z,t) dy = e i ^U^( t ^+^+ 5 ^^-^ TT 5 ( Xj - Zj ) , (2.34) 



a 



[ H(x,y,t)G(y,z,t) dy = e-^ifriWto+^iWX**-**) ff 5 ( x . _ z .) 
V fJi 

where 5 (x) is the Dirac delta function with respect to the space coordinates. 

Next, we have 

Uh 1 (t)U H (t)^j(x,0) = U H 1 (t)VM) 
H(x,y,t) ip(y,t) dy 



(2.35) 



H(x,y,t) I J ^G H (y,z,t) ij)(z,0) dzj dy 
H (x,y,t)G H (y,z,t) dy) ip(z,0) dz 
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» d 

= / e-E^t^Wfe+^+^flK^r^) TT $ ( Xj - Zj ) $ (z, 0) dz 
= ip(x,0) , 

or U H l {t) Uh [t) = I. A formal proof of the second relation Uh (t) U H X (t) = I is similar. The rest 
of the statement follows for a standard duality argument, see jH]- 

3. This claim is a consequence of (1) multiplying the equation 

c (t) 

iifrt = -a(t)ip xx + b(t)x 2 ip - f(t)xip - i(c(t)x - g(t))d x ip - i—^-^ 

by ip (2) integrating in the space variable, and (3) equating the imaginary parts of both sides to 
obtain 

Re J Mdx=~J {c{t)x-g{t))d x {\4>\ 2 )dx- l -J c(t) \^\ 2 dx; (2.36) 

the estimate now follows from the solution being smooth. 

We introduce the integral operator Uh {t, s) = Uh (t) U H l (s) by 



Uh (t) U^ (s) i> (x, s) = / G (x, y, t, s) $ (y, s) dy (2.37) 
with the kernel given by 

G(x,y,t,s)— I Gh {x, z,t) H (z,y, s) dz. (2.38) 

Here, 



( d 1 Y /2 

G (x, y, t,s)=[ (2.39) 

d 

x yj exp (i [ocj (t) x 2 - ctj (s) y 2 + Sj (t) Xj - Sj (s) yj + Kj (t) - Kj (s))) 

3=1 

/ (Pj (t) - 0, (s) y 3 + Sj (t) - Sj (s)) 2 



x exp 

3=1 



Ai ( 7j {t) - 7j (s)) 



\U H (t, s)ip(x,s)\ 



< 



n 



G(x,y,t,8)ip(y,s) dy 
1 



d/2 



i V 471 ^ (*) ^ ( s ) Wi ( s ) - 7j (*)) 



|^(y,s)| dy. 



Thus, the estimate (12.121) holds. 

4. The uniqueness of the solution and continuous dependence on the initial data follows by 
standard arguments using estimates of the type of (12. lip . 
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In this section we will discuss Strichartz type estimates for the operator Uh- 

Definition 1. We say that the exponent pair (q,r) is a — admissible ifq,r > 2, (q,r,o~) ^ (2, oo, 1) 
and 

I a a 

- + -<«■ 
q r 2 

If equality holds we say that (q, r) is sharp a -admissible, otherwise, we say that (q, r) is nonsharp 
a -admissible. Note, in particular, that when a > 1 the endpoint 



a 



is sharp a-admissible. 



The following inequalities are known as Strichartz estimates studied by Strichartz, Ginibre, Velo, 
Keel and Tao and others, see [31] for the following version. 

Lemma 4. IfU(t) ://"—>■ L 2 (x), where H is a Hilbert space and 

• For all t > and f G H, we have 

\\U(t)f\\ Ll <\\f\\ H , 

• For some a > 0, t ^ s and g G ^(x), 



\Mt)u( s y g \ 



oo<\t~ S\ 



Then 



\Mt)f\\rt 
(U(s))* F(s)ds 

U(t) (U(s))* F(s)ds 



L 2 



<t 



< 



< 



\F\ 
\F\ 



H 



'1 TV 1 



(3.1) 

(3.2) 

(3.3) 
(3.4) 

(3.5) 



hold for all sharp a-admissible exponent pairs (q,r), (q,r). 

In order to obtain inequalities (13. 3p . (13. 4p . (13. 5p for our operator Uh, we observe that the inequal- 
ities are valid if we replace (13. 2 p by 

\\U(t)(U(s)yF(s)\\ L ~<w(t-sy\\F(s)\\ Li , weLl (3.6) 

We can just mimic the proof of Lemma 3, see [56] and the adaptation in [5] done for the case Uy. 
For our case we deal with the inequality (13.61) and avoid the use of semigroup properties that we 
do not have. More specifically we obtain: 

Lemma 5. If Unit) satisfies k3. 6\) . then for any T G M +7 

\\U H (t)f\y tL r < c q \\w-l ( _ 2T , 2T) \\2 2 \\f\\ Lh (3.7) 



(U H (s))* F(s)ds 



L 2 



< Cg\\w-1^2T,2T)\r L 2 2\\F 



(3. 
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U H (t) (U H (s))* F(s)ds 



s<t 



< C qt q\\lV-l 



173 



L q t {{-T,T);Lr x ) 



L(-2T,2T)|| £ 2||-r \ \ L f L r>, 



(3.9) 



holding for all sharp a-admisible exponent pairs (q,r), (q,r), where ji, 72 ,73 are constants. 
For the sake of clarity we outline the proof: 

First, observe that (13.41) implies (13. 3p by duality. Therefore we prove ()3l 



U H (s)*F(s)ds 



< \ \F\ 



Li' (R, L T ' 



L 2 



By the TT* method this can be implied by the following inequality: 



< ILFI 




< U H (s)*F(s),U H (tyG{t)dsdt > 
By symmetry, it suffices to prove 

\T(F n\\ < \ \F 



A\G\\, 



q' j T l ■ 



\t{f } G)\<\\f\\^ t ^\\g\\ l ^ 



(3.10) 



where 



\T(F,G)\ 




< UH(s)*F(s),U H {t)*G{t) > dsdt. 



s<t 



Since Uh is unitary (U^ = Uy- 1 ) on L 2 , by Holder's inequality and the energy estimate, we get 

\<u H ( s yF(s),u H (tyG(t)>\ < \\u H ( s yF(s)\\ Ll \\u H (tyG(t)\\Li (3.11) 

= \\F( S )\\ Li \\G( S )\\ L2x . (3.12) 
By assumption (13.61) and from the above estimate, we get 

\<U H {s)*F(s),U H {tyG(t)>\ < \\U H (t)U H ( S yF(s)\\ L °o\\G{t)\\ L i (3.13) 

< u(s-ty\\F(s)\\ Li \\G(s)\\ Li . (3.14) 
The last equality follows by Holder's inequality and 7 is coming from (13. 6p . 
Now by interpolating with 

\<U H ( S yF(s),U H (tyG(t) >\ < \\F(s)\\ Ll \\G(s)\\ Lh 

we obtain that if r' is defined by l/r' = l — 6 + 9/2, < 9 < 1, then 1 < r' < 2. If we denote by r 
the dual exponent of r' then 1/r = 9/2, and choosing q such that (q, r) is 7-admissible, we will also 
get that (q',r f ) is sharp 7-admissible, and that 7(1 — 9) = 2/q. Therefore, 

|< U H ( S yF(s), U H (tyG(t) >\ < u(t - s)l\\F(s)\\ Lr ,\\G(s)\\ LrJ . 

Since (g, r) is sharp 7-admissible (and it is not an endpoint), then we can apply weak young 
inequality, obtaining 

J j * \<U H {s)*F{s),U H {t)*G{t) >\dsdt < J J u{t - s)-«\\F{s)\\ LT j\\F{s)\\ G r,dsdt 

< \Wt-s)\\ 



\F\\ T <i' jA\G\ 



where 2/q + l/q' + l/q' = 2, and since u is in weak L , (13 Ah follows. 
Now to prove (13.51) we can proceed as in [56] (Section 7). 





= -1, 




= 0, 




= +1, 
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The following result gives us conditions on operator Uh to obtain global in time Strichartz 
estimates: 

Theorem 2. 1. Consider the following restriction on the coefficients 

j>lt ffj . e{ _i,o,i } . (3.i5) 

T/ie evolution operator Uh associated to the Cauchy problem U.1\) - /[T72)) is given by 112.1 0\) where 

{sinh(oj jt) 
t, i/o-i = 0, (3.16) 
s\n(u> jt) 

and it satisfies 

( d i \ dn 

II^ML<(n^^j Ml.- ("7) 

Furthermore, if Sj = —1 /or some j, we have global in time Strichartz estimates. 

Proof. It is easy to see that (I3.17P follows from (12.121) and (13.161) . The global in time Strichartz 
estimates follow from Lemma 5 and observing that (as was pointed out in [5] (Section 2)) if Sf. = —1 
for some k and 5j = 1 for i ^ k (the worst of the possible cases) , then 

= c (w\ lm> + (^'n 2^4^) ' < 3 - i8 » 

is in (R), and from Lemma 5 we obtain global in time Strichartz estimates. □ 



4. DAMPED HARMONIC NONLINEAR SCHRODINGER EQUATION 

All linear exactly solvable models which we discussed in Sections 1 and 2 are of interest in a 
general treatment of the nonlinear time-dependent Schrodinger equation, see [26] , [28] . [40] . |46j. 
In this Section we apply the results of the last sections to the study of the nonlinear version of 
equation (ll.ip - (11.2p (with algebraic nonlinearity) : 

> d i = 4 A " + E (^ 2 " - W)v + i9i (t)^ - (^g - «)) + hwr'u (4.1) 

u(x,0) = u (x). (4.2) 

It includes the damped harmonic nonlinear Schrodinger equation (11. 4p and the following well- 
known cases: 

• The nonlinear Schrodinger equation with zero potential 

f)ni j 

ij- = --Au + hu\u\ p - 1 . (4.3) 

L/C —i 
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The nonlinear Schrodinger equation with the quadratic potential possibly depending on time 

-- Au + J2^-x 2 j u + hu\u\ p - 1 . (4.4) 



du 

% 



dt 2 ^2 

3=1 



• The Gross-Pit aevskii equation 

du \ TfL 

i— = --Am + —u(t) 2 x 2 u - f{t)xu + hu\u\ p ~ l , (4.5) 
see [60], [39]. 

These equations have different applications, for example, the propagation of waves, optical trans- 
mission lines with online modulators, propagation of light beams in nonlinear media with a gradient 
of the refraction index, or in the theory of Bose-Einstein condensate in trapped gases [16] , [60] . The 
following local well-posedness result in L 2 x (R d ) in the subcritical sense is a consequence of Strichartz 
estimates and a boostrap argument |56j : 

Proposition 1. Let p be an L 2 — subcritical exponent (0 < p — 1 < 4/d), h = ±1. Then for 
any R > there exists T > such that for all u E L 2 x (R d ) in the ball B R = {u E Ll(R d ) : 
||wo|lL2( R d) < R} there exists a unique strong L 2 , solution u to ( |^.i| )-( [7^P in the space S°([—T,T] x 
R d ) cC?Ll([-T,T\ x R d ). 

Sketch of the proof: Since we obtained Strichartz estimates, see Section 3, we can follow the proof 
for the case of the nonlinear Schrodinger equation without potential, see [56], [8] or [51]. 

Remark 4. (Time dependent nonlinearity) One can also consider a generalization of the nonau- 
tonomous Schrodinger equation 

= —AU+J2 ( bj ^*> - fMxju + i 9j (t)^- - (2x^ - u\ ) + h(t)\ur 1 u (4.6) 



u(x,0) = u (x), (4.7) 
considering for example h E C°°(R; R). This equation will include the nonautonomous Schrodinger 
equation [S] 

r)ii 1 

i^- = -=-Au + h{t)u\u\ p - 1 . (4.8) 
at 2 
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